Machine-learning analysis reveals an important role for negative selection in shaping cancer aneuploidy landscapes.

BACKGROUND: Aneuploidy, an abnormal number of chromosomes within a cell, is a hallmark of cancer. Patterns of aneuploidy differ across cancers, yet are similar in cancers affecting closely related tissues. The selection pressures underlying aneuploidy patterns are not fully understood, hindering our understanding of cancer development and progression. RESULTS: Here, we apply interpretable machine learning methods to study tissue-selective aneuploidy patterns. We define 20 types of features corresponding to genomic attributes of chromosome-arms, normal tissues, primary tumors, and cancer cell lines (CCLs), and use them to model gains and losses of chromosome arms in 24 cancer types. To reveal the factors that shape the tissue-specific cancer aneuploidy landscapes, we interpret the machine learning models by estimating the relative contribution of each feature to the models. While confirming known drivers of positive selection, our quantitative analysis highlights the importance of negative selection for shaping aneuploidy landscapes. This is exemplified by tumor suppressor gene density being a better predictor of gain patterns than oncogene density, and vice versa for loss patterns. We also identify the importance of tissue-selective features and demonstrate them experimentally, revealing KLF5 as an important driver for chr13q gain in colon cancer. Further supporting an important role for negative selection in shaping the aneuploidy landscapes, we find compensation by paralogs to be among the top predictors of chromosome arm loss prevalence and demonstrate this relationship for one paralog interaction. Similar factors shape aneuploidy patterns in human CCLs, demonstrating their relevance for aneuploidy research. CONCLUSIONS: Our quantitative, interpretable machine learning models improve the understanding of the genomic properties that shape cancer aneuploidy landscapes.

Genetic underpinnings explored: OPA1 deletion and complex phenotypes on chromosome 3q29.

BACKGROUND: Copy number variations (CNVs) have emerged as significant contributors to the elusive genetic causality of inherited eye diseases. In this study, we describe a case with optic atrophy and a brain aneurysm, in which a de novo CNV 3q29 deletion was identified. CASE PRESENTATION: A 40-year-old female patient was referred to our department after undergoing aneurysm transcatheter arterial embolization for a brain aneurysm. She had no history of systemic diseases, except for unsatisfactory best-corrected visual acuity (BCVA) since elementary school. Electrophysiological tests confirmed the findings in retinal images, indicating optic nerve atrophy. Chromosomal microarray analysis revealed a de novo deletion spanning 960 kb on chromosome 3q29, encompassing OPA1 and six neighboring genes. Unlike previously reported deletions in this region associated with optic atrophy, neuropsychiatric disorders, and obesity, this patient displayed a unique combination of optic atrophy and a brain aneurysm. However, there is no causal relationship between the brain aneurysm and the CNV. CONCLUSION: In conclusion, the optic atrophy is conclusively attributed to the OPA1 deletion, and the aneurysm could be a coincidental association. The report emphasizes the likelihood of underestimating OPA1 deletions due to sequencing technology limitations. Recognizing these constraints, healthcare professionals must acknowledge these limitations and consistently search for OPA1 variants/deletions in Autosomal Dominant Optic Atrophy (ADOA) patients with negative sequencing results. This strategic approach ensures a more comprehensive exploration of copy-number variations, ultimately enhancing diagnostic precision in the field of genetic disorders.

An Optimized Genotyping Workflow for Identifying Highly SCRaMbLEd Synthetic Yeasts.

Synthetic Sc2.0 yeast strains contain hundreds to thousands of loxPsym recombination sites that allow restructuring of the Saccharomyces cerevisiae genome by SCRaMbLE. Thus, a highly diverse yeast population can arise from a single genotype. The selection of genetically diverse candidates with rearranged synthetic chromosomes for downstream analysis requires an efficient and straightforward workflow. Here we present loxTags, a set of qPCR primers for genotyping across loxPsym sites to detect not only deletions but also inversions and translocations after SCRaMbLE. To cope with the large number of amplicons, we generated qTagGer, a qPCR genotyping primer prediction tool. Using loxTag-based genotyping and long-read sequencing, we show that light-inducible Cre recombinase L-SCRaMbLE can efficiently generate diverse recombination events when applied to Sc2.0 strains containing a linear or a circular version of synthetic chromosome III.

A chromosome-level genome assembly for the amphibious plant Rorippa aquatica reveals its allotetraploid origin and mechanisms of heterophylly upon submergence.

The ability to respond to varying environments is crucial for sessile organisms such as plants. The amphibious plant Rorippa aquatica exhibits a striking type of phenotypic plasticity known as heterophylly, a phenomenon in which leaf form is altered in response to environmental factors. However, the underlying molecular mechanisms of heterophylly are yet to be fully understood. To uncover the genetic basis and analyze the evolutionary processes driving heterophylly in R. aquatica, we assembled the chromosome-level genome of the species. Comparative chromosome painting and chromosomal genomics revealed that allopolyploidization and subsequent post-polyploid descending dysploidy occurred during the speciation of R. aquatica. Based on the obtained genomic data, the transcriptome analyses revealed that ethylene signaling plays a central role in regulating heterophylly under submerged conditions, with blue light signaling acting as an attenuator of ethylene signal. The assembled R. aquatica reference genome provides insights into the molecular mechanisms and evolution of heterophylly.

Chromosome-level genome assembly of the yellow-cheek carp Elopichthys bambusa.

Yellow-cheek carp (Elopichthys bambusa) is a typical large and ferocious carnivorous fish endemic to East Asia, with high growth rate, nutritional value and economic value. In this study, a chromosome-level genome of yellow-cheek carp was generated by combining PacBio reads, Illumina reads and Hi-C data. The genome size is 827.63 Mb with a scaffold N50 size of 33.65 Mb, and 99.51% (823.61 Mb) of the assembled sequences were anchored to 24 pseudo-chromosomes. The genome is predicted to contain 24,153 protein-coding genes, with 95.54% having functional annotations. Repeat elements account for approximately 55.17% of the genomic landscape. The completeness of yellow-cheek carp genome assembly is highlighted by a BUSCO score of 98.4%. This genome will help us understand the genetic diversity of yellow-cheek carp and facilitate its conservation planning.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers.

BACKGROUND: Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei. RESULTS: Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions. CONCLUSIONS: Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

Chromosome-level genome assembly of marine diatom Skeletonema tropicum.

Skeletonema tropicum is a marine diatom of the genus Skeletonema that also includes many well-known species including S. marinoi. S. tropicum is a high temperature preferring species thriving in tropical ocean regions or temperate ocean regions during summer-autumn. However, mechanisms of ecological adaptation of S. tropicum remain poorly understood due partially to the lack of a high-quality whole genome assembly. Here, we report the first high-quality chromosome-scale genome assembly for S. tropicum, using cutting-edge technologies including PacBio single molecular sequencing and high-throughput chromatin conformation capture. The assembled genome has a size of 78.78 Mb with a scaffold N50 of 3.17 Mb, anchored to 23 pseudo-chromosomes. In total, 20,613 protein-coding genes were predicted, of which 17,757 (86.14%) genes were functionally annotated. Collinearity analysis of the genomes of S. tropicum and S. marinoi revealed that these two genomes were highly homologous. This chromosome-level genome assembly of S. tropicum provides a valuable genomic platform for comparative analysis of mechanisms of ecological adaption.

Chromosome Transplantation: Opportunities and Limitations.

There are thousands of rare genetic diseases that could be treated with classical gene therapy strategies such as the addition of the defective gene via viral or non-viral delivery or by direct gene editing. However, several genetic defects are too complex for these approaches. These "genomic mutations" include aneuploidies, intra and inter chromosomal rearrangements, large deletions, or inversion and copy number variations. Chromosome transplantation (CT) refers to the precise substitution of an endogenous chromosome with an exogenous one. By the addition of an exogenous chromosome and the concomitant elimination of the endogenous one, every genetic defect, irrespective of its nature, could be resolved. In the current review, we analyze the state of the art of this technique and discuss its possible application to human pathology. CT might not be limited to the treatment of human diseases. By working on sex chromosomes, we showed that female cells can be obtained from male cells, since chromosome-transplanted cells can lose either sex chromosome, giving rise to 46,XY or 46,XX diploid cells, a modification that could be exploited to obtain female gametes from male cells. Moreover, CT could be used in veterinary biology, since entire chromosomes containing an advantageous locus could be transferred to animals of zootechnical interest without altering their specific genetic background and the need for long and complex interbreeding. CT could also be useful to rescue extinct species if only male cells were available. Finally, the generation of "synthetic" cells could be achieved by repeated CT into a recipient cell. CT is an additional tool for genetic modification of mammalian cells.

Chromosome-scale assembly of the streamlined picoeukaryote Picochlorum sp. SENEW3 genome reveals Rabl-like chromatin structure and potential for C4 photosynthesis.

Genome sequencing and assembly of the photosynthetic picoeukaryotic Picochlorum sp. SENEW3 revealed a compact genome with a reduced gene set, few repetitive sequences, and an organized Rabl-like chromatin structure. Hi-C chromosome conformation capture revealed evidence of possible chromosomal translocations, as well as putative centromere locations. Maintenance of a relatively few selenoproteins, as compared to similarly sized marine picoprasinophytes Mamiellales, and broad halotolerance compared to others in Trebouxiophyceae, suggests evolutionary adaptation to variable salinity environments. Such adaptation may have driven size and genome minimization and have been enabled by the retention of a high number of membrane transporters. Identification of required pathway genes for both CAM and C4 photosynthetic carbon fixation, known to exist in the marine mamiellale pico-prasinophytes and seaweed Ulva, but few other chlorophyte species, further highlights the unique adaptations of this robust alga. This high-quality assembly provides a significant advance in the resources available for genomic investigations of this and other photosynthetic picoeukaryotes.

A historical sequence deletion in a commonly used Bacillus subtilis chromosome integration vector generates undetected loss-of-function mutations.

Since the 1980s, chromosome-integration vectors have been used as a core method of engineering Bacillus subtilis. One of the most frequently used vector backbones contains chromosomally derived regions that direct homologous recombination into the amyE locus. Here, we report a gap in the homology region inherited from the original amyE integration vector, leading to erroneous recombination in a subset of transformants and a loss-of-function mutation in the downstream gene. Internal to the homology arm that spans the 3' portion of amyE and the downstream gene ldh, an unintentional 227 bp deletion generates two crossover events. The major event yields the intended genotype, but the minor event, occurring in ~10 % of colonies, results in a truncation of ldh, which encodes lactate dehydrogenase. Although both types of colonies test positive for amyE disruption by starch plating, the potential defect in fermentative metabolism may be left undetected and confound the results of subsequent experiments.

Characterization and Mapping of a Rolling Leaf Mutant Allele rlT73 on Chromosome 1BL of Wheat.

Leaf rolling is regarded as an important morphological trait in wheat breeding. Moderate leaf rolling is helpful to keep leaves upright and improve the photosynthesis of plants, leading to increased yield. However, studies on the identification of genomic regions/genes associated with rolling leaf have been reported less frequently in wheat. In this study, a rolling leaf mutant, T73, which has paired spikelets, dwarfism, and delayed heading traits, was obtained from a common wheat landrace through ethyl methanesulfonate mutagenesis. The rlT73 mutation caused an increase in the number of epidermal cells on the abaxial side and the shrinkage of bulliform cells on the adaxial side, leading to an adaxially rolling leaf phenotype. Genetic analysis showed that the rolling leaf phenotype was controlled by a single recessive gene. Further Wheat55K single nucleotide polymorphism array-based bulked segregant analysis and molecular marker mapping delimited rlT73 to a physical interval of 300.29-318.33 Mb on the chromosome arm 1BL in the Chinese Spring genome. We show that a point mutation at the miRNA165/166 binding site of the HD zipper class III transcription factor on 1BL altered its transcriptional level, which may be responsible for the rolling leaf phenotype. Our results suggest the important role of rlT73 in regulating wheat leaf development and the potential of miRNA-based gene regulation for crop trait improvement.

Chromosome-level genome assembly of Hippophae tibetana provides insights into high-altitude adaptation and flavonoid biosynthesis.

BACKGROUND: As an endemic shrub of the Qinghai-Tibetan Plateau (QTP), the distribution of Hippophae tibetana Schlecht. ranges between 2800 and 5200 m above sea level. As the most basal branch of the Hippophae genus, H. tibetana has an extensive evolutionary history. The H. tibetana is a valuable tree for studying the ecological evolution of species under extreme conditions. RESULTS: Here, we generated a high-quality chromosome-level genome of H. tibetana. The total size of the assembly genome is 917 Mb. The phylogenomic analysis of 1064 single-copy genes showed a divergence between 3.4 and 12.8 Mya for H. tibetana. Multiple gene families associated with DNA repair and disease resistance were significantly expanded in H. tibetana. We also identified many genes related to DNA repair with signs of positive selection. These results showed expansion and positive selection likely play important roles in H. tibetana's adaptation to comprehensive extreme environments in the QTP. A comprehensive genomic and transcriptomic analysis identified 49 genes involved in the flavonoid biosynthesis pathway in H. tibetana. We generated transgenic sea buckthorn hairy root producing high levels of flavonoid. CONCLUSIONS: Taken together, this H. tibetana high-quality genome provides insights into the plant adaptation mechanisms of plant under extreme environments and lay foundation for the functional genomic research and molecular breeding of H. tibetana.

Chromosome-level genome provides insights into environmental adaptability and innate immunity in the common dolphin (delphinus delphis).

The common dolphin (Delphinus delphis) is widely distributed worldwide and well adapted to various habitats. Animal genomes store clues about their pasts, and can reveal the genes underlying their evolutionary success. Here, we report the first high-quality chromosome-level genome of D. delphis. The assembled genome size was 2.56 Gb with a contig N50 of 63.85 Mb. Phylogenetically, D. delphis was close to Tursiops truncatus and T. aduncus. The genome of D. delphis exhibited 428 expanded and 1,885 contracted gene families, and 120 genes were identified as positively selected. The expansion of the HSP70 gene family suggested that D. delphis has a powerful system for buffering stress, which might be associated with its broad adaptability, longevity, and detoxification capacity. The expanded IFN-α and IFN-ω gene families, as well as the positively selected genes encoding tripartite motif-containing protein 25, peptidyl-prolyl cis-trans isomerase NIMA-interacting 1, and p38 MAP kinase, were all involved in pathways for antiviral, anti-inflammatory, and antineoplastic mechanisms. The genome data also revealed dramatic fluctuations in the effective population size during the Pleistocene. Overall, the high-quality genome assembly and annotation represent significant molecular resources for ecological and evolutionary studies of Delphinus and help support their sustainable treatment and conservation.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum.

BACKGROUND: Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism. RESULTS: We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening. CONCLUSION: Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein.

BACKGROUND: Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood. RESULTS: To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes. CONCLUSIONS: This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

The scales, mechanisms, and dynamics of the genome architecture.

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

Influence of the number of washings for embryos on non-invasive preimplantation chromosome screening results.

Objective: To explore the effect of varying numbers of embryo washings prior to blastocyst formation in non-invasive preimplantation chromosome screening (NICS) on the accuracy of NICS results. Methods: In this study, 68 blastocysts from preimplantation genetic testing (PGT)-assisted pregnancy were collected at our institution. On the fourth day of embryo culture, the embryos were transferred to a new medium for blastocyst culture and were washed either three times (NICS1 group) or ten times (NICS2 group). A trophectoderm (TE) biopsy was performed on the blastocysts, and the corresponding embryo culture media were collected for whole genome amplification (WGA) and high-throughput sequencing. Results: The success rate of WGA was 100% (TE biopsy), 76.7% (NICS1 group), and 89.5% (NICS2 group). The success rate of WGA in embryo medium on days 5 and 6 of culture was 75.0% (33/44) and 100% (24/24), respectively. Using TE as the gold standard, the karyotype concordance rate between the results of the NICS1 and NICS2 groups' embryo culture medium samples and TE results was 43.5% (10/23) and 73.5% (25/34), respectively. The sensitivity and specificity of detecting chromosomal abnormalities were higher in the NICS2 group than in the NICS1 group when TE was used (83.3% vs 60.0%; 62.5% vs 30.8%, respectively). The false-positive rate and false-negative rate (i.e., misdiagnosis rate and missed diagnosis rate, respectively) were lower in the NICS2 group than in the NICS1 group (37.5% vs 69.2%; 16.7% vs 40.0%, respectively). Conclusion: The NICS yielded favorable results after ten washings of the embryos. These findings provide a novel method for lowering the amount of cell-free DNA contamination from non-embryonic sources in the medium used for embryo development, optimizing the sampling procedure and improving the accuracy of the NICS test.

The chromosome-level genome assembly of the giant dobsonfly Acanthacorydalis orientalis (McLachlan, 1899).

Acanthacorydalis orientalis (McLachlan, 1899) (Megaloptera: Corydalidae) is an important freshwater-benthic invertebrate species that serves as an indicator for water-quality biomonitoring and is valuable for conservation from East Asia. Here, a high-quality reference genome for A. orientalis was constructed using Oxford Nanopore sequencing and High throughput Chromosome Conformation Capture (Hi-C) technology. The final genome size is 547.98 Mb, with the N50 values of contig and scaffold being 7.77 Mb and 50.53 Mb, respectively. The longest contig and scaffold are 20.57 Mb and 62.26 Mb in length, respectively. There are 99.75% contigs anchored onto 13 pseudo-chromosomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the completeness of the genome assembly is 99.01%. There are 10,977 protein-coding genes identified, of which 84.00% are functionally annotated. The genome contains 44.86% repeat sequences. This high-quality genome provides substantial data for future studies on population genetics, aquatic adaptation, and evolution of Megaloptera and other related insect groups.

Identification and functional analysis of a chromosome 2D fragment harboring TaFPF1 gene with the potential for yield improvement using a late heading wheat mutant.

KEY MESSAGE: A chromosome fragment influencing wheat heading and grain size was identified using mapping of m406 mutant. The study of TaFPF1 in this fragment provides more insights into wheat yield improvement. In recent years, wheat production has faced formidable challenges driven by rapid population growth and climate change, emphasizing the importance of improving specific agronomic traits such as heading date, spike length, and grain size. To identify potential genes for improving these traits, we screened a wheat EMS mutant library and identified a mutant, designated m406, which exhibited a significantly delayed heading date compared to the wild-type. Intriguingly, the mutant also displayed significantly longer spike and larger grain size. Genetic analysis revealed that a single recessive gene was responsible for the delayed heading. Surprisingly, a large 46.58 Mb deletion at the terminal region of chromosome arm 2DS in the mutant was identified through fine mapping and fluorescence in situ hybridization. Thus, the phenotypes of the mutant m406 are controlled by a group of linked genes. This deletion encompassed 917 annotated high-confidence genes, including the previously studied wheat genes Ppd1 and TaDA1, which could affect heading date and grain size. Multiple genes in this region probably contribute to the phenotypes of m406. We further investigated the function of TaFPF1 using gene editing. TaFPF1 knockout mutants showed delayed heading and increased grain size. Moreover, we identified the direct upstream gene of TaFPF1 and investigated its relationship with other important flowering genes. Our study not only identified more genes affecting heading and grain development within this deleted region but also highlighted the potential of combining these genes for improvement of wheat traits.